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NOMENCLATURE 


A area 

a speed of sound 

B jet layer parameter 

C v mass flow coefficient 

d diameter 

F(c) function describing resonant properties of the generator 
f(M) function describing acoustical properties of nozzle 
f frequency 

fp diaphragm-reed friction coefficient 

g 0 acceleration of gravity, sea-level value 

H geopotential altitude 

h 0 ram-air jet width 

I flow inertia parameter 

i = AT 

K non-dimensior*al cavity volume 

kp diaphragm-reed spring constant 

A distance from ram-air jet to tube 

L tube length 

M( ) function describing resonant characteristics of the tube 
M 0 molecular weight of air 

Mi generator flight Mach number 

mp diaphragm-reed effective mass 

m mass flow 

p pressure 

q ram-air jet velocity 

R reflection coefficient 



Nomenclature, (Continued) 


R* universal gas constant 

Rj radius of curvature, jet layer 

liid 

R e * — , Reynolds number 

An 

r 0 effective earth radius 

T temperature 

t time 

U velocity 

V volume 

v velocity component across the flow direction 

W jet layer perimeter 

x downstream distance 

xq diaphragm-reed position 

1 geometrical altitude 

z = ( 0+1 to) complex frequency 

a,p non-dimensional parameters 

V specific heat ratio 

C non-dimensional complex frequency 

t) diaphragm-reed non-dimensional parameter 

p viscosity 

itij 

W ” “ 

my 

m H 

PH = — 

my 

a real part complex frequency 



Nomenclature, (Continued) 


$,q diaphragm-reed non-dimensional parameters 

m Imaginary part complex frequency 

SUBSCRIPTS AND SUPERSCRIPTS 

( ) 0 corresponds to stagnation conditions 

( )\ conditions In front of normal shock 

( )g conditions behind the normal shock 

(') rate of flow 

( )' perturbation quantity 

( )„ free-stream flow conditions 

( ) c corresponds to cavity conditions (See Figure 20) 

{ )d relates to diaphragm-reed system 

( )e relates to conditions at the end of the tube 

( )n relates to flow conditions across exit holes (Figure 20) 

( )j corresponds to annular ram-air jet 

( )s relates to static flow conditions at the generator exit holes 
( )y relates to resonance tube conditions 

( ) w working equilibrium flow conditions 




ABSTRACT 


A combined experimental and analytical investigation was performed to 
understand the various fluid processes associated with the conversion of flow 
energy into electric power In a fluidic generator. Experiments were performed 
under flight-simulated laboratory conditions and results were compared with 
those obtained in the free-flight conditions. From these measurements, It was 
concluded that the mean mass flow critically controlled the output of the 
fluidic generator. Cross-correlation of the outputs of transducer data 
Indicated the presence of a standing wave In the tube; therefore the mechanism 
of oscillation was an acoustic resonance tube phenomenon and did not represent 
the Helmholtz type of oscillation phenomenon. 

On the basis of experimental results, a linearized model was constructed 
coupling the flow behavior of the jet, the jet-layer, the tube, the cavity, 
and the holes of the fluidic generator. The analytical results also showed 
that the mode of the fluidic power generator is an acoustical resonance 
phenomenon with the frequency of operation given by f« a/4L, where f is the 
frequency of jet swallowing, a is the average speed of sound in the tube, and 
L is the length of the tube. Analytical results further Indicated that 
oscillations in the fluidic generator are always damped and consequently there 
Is a forcing of the system in operation. 


x 



I. INTRODUCTION 


Periodic oscillations In jet flows Impinging on solid surfaces have been 
observed over a range of flow conditions and geometrical variations (Reference 
1). The oscillation phenomenon In such jet flows is very similar to that 
observed In axlsymmetric cavity flows. Both experimental and anu\vt1c&T 
investigations of flows over cavities by V. Sarohla (References 2 to 6) have 
shown that the self-sustained oscillations are caused by the Instability of 
the mean velocity profile and that they are amplified along the shear layer 
flow over the <javity. These previous Investigations (Reference 5) further 
showed that the flow conditions downstream of the cavity Influence the 
oscillation phenomenon as well as the shift of flow oscillations from one mode 
to another. The particular mode of the oscillation is determined by the phase 
of the propagating disturbance which attains the maximum value of the 
Integrated amplification along the cavity shear layer. 

Experiments performed on axlsymmetric cavity flow oscillations showed the 
oscillation phenomenon depends on the external acoustic excitation (Reference 
7). The jet edge-tone configuration in the fluidic electric power generator, 
as sketched in Figure 1, is similar to one in the axisymmetric cavity flow. 
However, the jet edge-tone oscillations will be modified by the resonance 
tube, reed-diaphr«gm system, and cavity surrounding the fluidic generator. 

In flight, as indicated in Figure 1, air enters the fluidic power genera¬ 
tor through the port located in the nose and leaves through the exhaust ports 
spaced uni formally around the circumferences of the ogive. During the passage 
of the air through the generator periodic air flow oscillations are produced 
by the jet interaction with the resonance tube. These flow oscillations are 
transmitted through a diaphragm and rod to reed that switches magnetic flux 
within the coil, thereby inducing a voltage at the coil terminals. 


-1- 



Little Is known, however, about the manner in which the amplitude and the 
frequency of these flow oscillations Influence the conversion of flow energy 
Into electric power In the system shown In Figures 1 and 7. The understanding 
of various flow Interactions that occur In the fluidic power generator Is very 
essential to the efficiency of conversion of flow energy into electric power. 
The Investigation reported here was designed to advance the understanding of 
these unanswered fundamental questions. 

Laboratory measurements were made to relate to the free-fllght 
performance of the fluidic generator. Based on the various underlying fluid 
processes associated with the production of oscillations In the generator, an 
analytic model which can predict the operation and the influence of 
configuration changes on generator performance was constructed. This model is 
the subject of the present report. The Information obtained in this study has 
led to methods to Improve the overall efficiency of the conversion of flow 
energy to electric power generation. 



II. EXPERIMENTAL ARRANGEMENTS AND MEASUREMENTS 


To enhance our understanding of the various flow processes responsible 
for the conversion of flow energy into electric power, two experimental setups 
were constructed. The bench test setup is shown in Figure 2. This facility 
could supply compressed air at various stagnation temperatures up to 700°F. 

The air was heated by passing it through an electric heater with a temperature 
control Variac as shown in Figure 2. This controlled-temperature air was 
accelerated through a convergent nozzle which was flushed with the inlet of 
the fluidic power generator. The total temperature and pressure in the 
settling chamber, along with static pressure at the inlet were measured and 
utilized to compute the rate of mean mass flow entering the fluidic generator. 

The fluidic generator output, as indicated in Figure 3, was monitored on 
the oscilloscope as well as on the rms meter. An environment chamber was also 
built to vary the ambient static pressure around the fluidic generator. The 
.cavity pressure p c and ambient pressure p® were "hionitored. Figure 3 also 
indicates the various electric outputs of the fluidic generator, such as the 
regulator output, the generator output, and the frequency output, which were 
recorded throughout this investigation. , 

Since the bench test setup, as shown in Figure 2, did not simulate the 
air flow around the generator, experiments were performed by placing the 
generator in the free-jet experimental test setup shown in Figure 4. The 
total stagnation inlet pressure, cavity pressure p c , and the pressure at the 
exit holes, along with the output of the fluidic generator were measured at 
various flight-simulated flow conditions. 

To determine the mechanism of operation of the fluidic generator, three 
1/8-in. fast-response pressure transducers with a rise time of 2psec and 


-3- 




a frequency response from 2 to 40,000 H 2 were employed. One transducer was 
located at the base and the second was located flushed with the mouth of the 
tuhe to measure time varying pressure In the tube* The outputs of these 
transducers were cross-correlated on an all-digital-correlator. The third 
transducer was utilized to measure the cavity pressure fluctuations. Figures 
5 and 6 show the physical locations of these transducers In the fluidic 
generator. 

To determine the dynamic response and the natural frequency of operation 
of the generator diaphragm-reed system, an accelerometer was mounted to the 
reed as shown in Figure 3. The output of this accelerometer was recorded on 
an rms meter and displayed on an oscilloscope. 
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III. EXPERIMENTAL RESULTS 


3.1 Calculations of Laboratory Conditions to Simulate Free-Flight Conditions 
Given the altitude Z, the static conditions in front of the fluidic 
generator (denoted by subscript 1 in Figure 7) can be calculated using the 
U.S. Standard Atmosphere (Reference 8). To simplify calculations we use the 
geopotential altitude H, which equates the work required to lift a unit mass 
(under varying gravitational field) to Z with the work required to lift the 
same unit mass to H through a region with uniform gravity g 0 « 9.80655 m/s 3 . 


Transformation: 
From Z + H 


r 0 Z 

H * - ; r 0 « 6,356,766m 

r 0 +Z 


1 m. 15 - 0.006511 0<H<11,000m 

216.65 11,000<H<20,000 

216.65 + 0.001 (H-20,000) 20,000<H<30,000 

-5.2559 

0<H< U,Q00m 


288.15 

t 101,325.0 (-) 


( 1 ) 


22,632.1 exp [-1.5772 (11-11,000)10^] 

ll,000m<H<20,000m 

216.65 34.163 

[5,473.1 ( —-) 20,000m<H<30,00Qm 

T ra 


( 2 ) 


M 0 P* Po, 

P* ■ -T - ** 3.4837 r 

R* T ro T^ 

where M Q is the molecular weight and R* is the universal gas constant. In 

2 3 

these formulas T w is obtained in K, p w in newton per m , and p w in kg/m * 

,_ R* 

The sonic velocity at point (1) a® ^vvRV, for air y * 1.4, — * 287.05 

% 


joule/(kg K). 
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a„ a V^Ol.87 T„ - 20.046 Vt^, t„ l n kelvins 
U 

Mi * — , U Is the fluidic generator velocity In m/s. 


(3) 


Moving across the shock to point (2) 
p 2 2 y M* - (y-1) 7 -1 


PI 


Y+l 


P02 /(y +1 ) M 1 \ Y/Y / Y+l / m ^\ / 6 \ 

pT^ - \— - ) °\—) (—) 


6 

Y+l 


7 


2 Y Mi-( Y -1) 


2 (Y-l)Mj +2 M? + 5 


Z 

7Mi-l 


M2 


?. r Mi - (y-1) 7Mi -1 


(4) 


T 2 (2 Y «i - (y-1)][(y- 1) Hf +2] (7M*-1)(M^ + 5) 

T 1 2 2 2 
(y+1) Mi 36Mi 

The Important variable is the driving pressure Ap, which is the 
difference of the stagnation pressure p 0 2 at the front of the fluidic 
generator minus the cavity pressure around the resonator (Ap a p 0 2 - Pc)* To 
transform free flight to the bench test. It was assumed that Ap = p 0 2 - Pc * 
P 3 - p m (see Figure 3). Recause of difficulties in estimating p c in free 
flight, a similar idea was tried by equating (p 0 - p s ) or p 0 2 - p$ = P3 - P® 
(see Figures 3 and 7). 

Therefore, by knowing the altitude Z» and the fluidic generator velocity 
U, Mi - U/a ro (Z) at that altitude can be calculated from relations (1) to (3) 
above. Using shock relation (4), the stagnation pressure in front of the 
generator p 0 2 and therefore the laboratory conditions needed to simulate the 
flight can be computed. A computer program was developed, which is attached 
in Appendix A along with some sample calculations. 
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3.2 Comparison of Bench Results With Free-Jet Results 


To determine the Influence of flight flow around the fluidic generator, 
bench test results (Figure 2) were compared with those obtained in the free-jet 
facility (Figure 4). The driving pressure (p 0 - Pc) was kept constant. The 
pressure P 3 In Figure 3 is assumed equal to the stagnation pressure p 0 in front 
of the generator In the free-jet facility shown in Figure 2. Typical results 
of such a comparison are shown in Figure 8 . As is evident in these results, 
for a given differential pressure p 0 - p„, the bench test output was always 
greater than that obtained by the free-jet test results. It should be noted 
that for this comparison the same generator was utilized. A comparison with 
two other generators was also made, and similar results were obtained. These 
results indicated that the fluidic generator output did not depend upon the 
differential pressure p 0 - po,. 

Upon careful analysis of other experimental data, i.e., cavity pressure p c 
as a function of the driving pressure p 0 - p®, it was observed that for a given 
(Po - Poo)* the cavity pressure p c was always high in the free-jet case as 
compared with the bench tests. These results are shown in Figure 9. 

Therefore, in Figure 10, the fluidic generator output was plotted as a function 
of (Po " Pc) anr! the agreement between free-jet and bench test results was 
excellent. 

From this study, it is therefore concluded that the most important 
parameter that controls the fluidic generator output is not p 0 -Poo but rather 
Po -Pc» the differential pressure between the total stagnation pressure in 
front of the generator and the cavity pressure. 
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3.3 Comparison of Laboratory Results With Flight Data 

Apart from getting the physical insight needed to model the fluid 
electric power generator from the above experiments, an attempt was made to 
correlate the laboratory tests with those obtained from flight data. The 
flight data was provided to JPL by Marry Diamond Laboratories. From the 
knowledge of the flight velocity and altitude, the laboratory flow conditions 
were calculated as discussed above in Section 3.1 and Appendix A. 

The results of such a comparison are shown in Figure 11. Since p c was 
not available in the flight data, the static pressure behind the oblique shock 
was equated to p c . As is evident from Figure 11, the comparison with the 
bench test results is very poor. The discrepancy in such a comparison may 
have been caused by the assumption that cavity pressure p c in flight is equal 
to p ro around the generator close to the holes. Flight measurements of p c are 
needed to enhance our understanding to better correlate the flight data with 
laboratory test results. Furthermore, it should be noted that the absolute 
pressures in the bench facility and in flight are not identical. 
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3.4 Mass Flow Measurements 


To better understand the performance of the operation of the fluidic 
generator, detailed measurements of the mass flow through the fluidic 
generator were determined as a function of various operating pressures. The 
influence of the exit hole area on generator performance was also determined. 

The calculations utilized to determine the mass flow rate through the 
fluidic device are shown in Appendix B. Figure 12 indicates the mass flow 
rate normalized with cavity pressure and plotted as a function of pressure 
ratio Ps/Pc where p s is the static pressure at the holes, which under condi¬ 
tions of no external flow is the same as atmospheric pressure p*,. Tables 1 - 4 
give the measured data along with the calculated mass flow rate and computed 
nozzle flow coefficients. The mass flow rate, and consequently, generator 
output drop for a given pressure ratio p s /p c when the holes are closed 
partially. However, if the mass flow rate is normalized with stagnation 
pressure p 0 , then the results shown in Figure 12 collapse into a single curve 
as shown in Figure 13. Results of Figures 12 and 13 clearly indicate the 
importance of the cavity pressure p c in controlling the mass flow rate through 
the generator. Any changes in the geometry of the holes will influence p c 
which in turn modifies the mass flow rate and the output of the fluidic 
generator. 

To determine the influence of the ambient pressure p s = p® in the bench 
test on mass flow rate through the generator, the environment chamber 
shown in Figure 3 was employed. Figure 14 shows the typical results of such 
an investigation. For all the exit holes opened, the mass flow rate m was 
independent of the ambient pressure for a given stagnation pressure p 0 and p c . 
These results strongly suggest that the mass flow rate under subsonic flow 
conditions through the generator critically depends upon the differential 
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pressure (p 0 - p c ). 

The Influence of the Reynolds number calculated on the nozzle velocity 
coefficient Is shown In Figure 15. These results showed that the changes In 
the exit hole area did not significantly modify the ram-jet flow responsible 
for generating pressure fluctuations In the table. Tables 1-4 give the 
output of the nozzle velocity coefficient' for various operating conditions* 

3.5 Mechanism of Fluidic Generator Operation 

To correctly model the operation of the fluidic generator, important 
information was needed to determine the mechanism responsible for pressure 
fluctuations in the tube. As shown in Figures 5 and 6, the pressure 
transducer output was cross-correlated to determine any standing waves in the 
tube. 

The instant pressure transducer output at the mouth and base of the 
resonance tube of the fluidic generator is shown in Figure 16, To mount the 
transducer at the base, the diaphragm was replaced by a solid plate. Figure 
16 also shows the auto and cross-correlations of these microphone signals. A$ 
is evident, the peak of the cross-correlation signal was delayed by 1/4 the 
time period of the pressure fluctuation frequency in the tube. This frequency 
is the fundamental tube resonance frequency given by f=a/4L, where L is the 
length of the tube. From these measurements, it is concluded that the 
pressure fluctuations, and consequently the oscillations of the diaphragm-reed 
system of the fluidic generator, result from an acoustic resonance phenomenon 
in the tube. These measurements did not support the Helmholtz type of 
oscillation phenomenon in the generator as has been commonly assumed in the 
past. 

The influence of mass flow rate on the output of the pressure 
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oscillations at the base of the resonance tube (so called since the frequency 
of pressure oscillations is at the fundamental resonance tube frequency) was 
studied. Results in Figure 17 indicated a linear relationship between the 
mass flow rate and the rms pressure fluctuations at the base of the cavity. 
This linear relationship suggests that the annular ram-air jet Is periodically 
swallowed by the tube at the frequency of resonance tube oscillation. This 
phenomenon, discussed in detail in Reference 9, Is called the regurgitant mode 
of resonance tube operation. Furthermore, since mass flow rate Is 
proportional to (p 0 -p c ) as discussed in Section 3.4, the magnitude of the 
pressure fluctuations at the base of the resonance tube is directly related to 
this differential pressure. 

3.6 Accelerometer Output 

Other experimental information needed to properly model the fluidic 
generator required us to determine the transient response of the fluidic 
diaphragm-reed subsystem of the fluidic generator. The accelerometer output 
was compared with the generator output at different pressures, as shown in 
Figure 18. Both outputs were at the resonance tube frequency. The transient 
response of the diaphragm-reed system is shown in Figure 19. The present 
measurements show the natural frequency of oscillation of the diaphragm-reed 
subsystem to be about 1412 Hz. The information shown in Figures 18 and 19 was 
utilized in an analytical model employed to compute the performance of this 
fluidic generator. 
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IV. ANALYTICAL RESULTS 


4.1 Governing Equations of Fluidic Generator Operation 

Since the operation of the fluidic generator performance critically 
depends upon the mass flow rate through the system, the governing equations 
which couple the annular ram-air jet, ring-tone jet-layer, cavity and 
resonance tube have been written in terms of nomenclature and Figure 20. The 
governing equations can be written as 

dl V /dp c \ . 

Ij — +1 I m* = - p' Annular Jet 

dt \dmj/ J c 


rimy' 

IT — + By m' = p' - p' Jet Layer 

dt T T c 


a M . 

— m' = - p' 
Ay T t 


Resonance Tube 


dp c ' a^ a? 

— = — m* = — (m' - m' + m') Cavity 

dt V c c V c J H T 


IH 


/ d Pc\ 

dt V dm / 

H w 


p ; 


Holes 


dp c 

where ( )' denotes acoustical oscillations and (—) denotes 

dmj yy 

the cavity pressure with the jet mass flow at the working (w) 
fluidic generator. For detailed derivation of the equations. 


the variation of 
point of the 
see Appendix C. 
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For solutions of the form e**, z 




/PT' 

“c 


K 

m' 

J 

■ 


M 






<y+1o>, In vector notation we write 


exp (zt) 


where 


z - a + 1(j t M 


2zL 

1 + R exp [-] 

a 


1 


- R exp [- 


2zL 

—3 

a 


One finds In matrix notation 



w 


0-1 0 

i"' • 

\l 0 0 


■ - \ 


h 

0 0 


Pc 

/d p c \ 


A 

Ih z H ~ ) 0 

Wh/ 


"j 

0 I T z + B T , 

a M 


"H / 

0 ~Kf j 


1 At / 


M 


Introducing nondlmensional variables 


where 


at 

Ax /dp c \ 

a j - ij f i 

>PjB irW 

T AT 

«H B p - 

A T /dPc\ 

» PH = -vzv 

a 'dmn 

at 

At 

«t * it - i 

■ (*t ■ n t 


z L V c itih 

C s — = — . pj = — ; PH “ — 

a AtL * ? 

tny 

L * Length of the resonance tube 
Ay = Area of the resonance tube 
a * Speed of sound = average 
The governing equation can br written 


10 KC -1 1 -1 \ 

/ ° T \ 

0 1 ocjC+Pj 0 0 ] 


f o c 

0 -1 0 «hC+Ph 0 


PJ 

-1 1 0 0 axC+PT, 


i w/ 

11 0 0 0 M / 

\ i/ 


For solution, the determinant should be zero 
0 KC -1 1 -1 

0 1 axC+Pj 0 0 

0-1 0 anC+PH 0 * o 

-1111 a T C+PX 
1 0 0 0 M 
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which results In the equation 

F(C)s M + atC+PT [K?(«jC+Pj)(aHC+PH) + («jC+Pj)+(aHC+PH)3 
+ (ajC+pj)(aHC+pH) • 0 (5) 


1-t+nC 


with M > 


(■*S) 


R * • 


l+*+nC 


,2 

paA 0 


(* 5 ) 


"*D a L /kj) 

<? “- » n *-. Q * - / — 

AjfD fp L a V mo 

where D denotes the diaphragm-reed conditions. 

For the solution of the above equation, we need to calculate the 
constants first. 

4.2 Calculations of Constants 

For these calculations, the dimensions of the fluidic generator have been 
used as follows: 

Diaphragm-Reed 

For these calculations, data In Figure 19 has been utilized, 
mo a 


fD L 


22.5 


C = 0.946 
<t> is unknown 
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Jet 


h Ay 

Aj L 

.5243 


pc 

.006 — 
m 

7tl.(5 

\P< 

•T 

/P< 
5-6 (- 
Vpr 

: 2/7 


Holes 


«K 


- for p c /p 0 < 0.528 

d H a T 


.0936 


A H L 


,2/7, 


pc 7 [!“(Ps/Pc) 3 Pc 
PH * t” — .— ■ ■ ■ v.i for — 


2-(P s /p c ) 


2/7 


Pc p£ 

Ph = 0.0006 — for — < 0.528 
m p c 


Jet Layer 


hAy 

a T = - * 0.0132 

AWL 


Po 


PT = 


P a 


a2 W2 


« 349 rn 


Cavity 


K - — = 14.5831 
A T L 

Independent Variables 

C * (a+1w) 1 = 1.095 x 10’ 4 (cr+lw) 


-16 


> 0.528 


> 0.528 



4.3 Numerical Results 


A computer program was developed to solve the equation (5) as developed 
In 3.1 above with the above physical constant. The computer program Is listed 
In Appendix D. The Input variables were the working conditions as follows: 
p 0 stagnation pressure in front of the oscillations 
p c cavity pressure 

p s static pressure close to the exit holes 
m mass flow rate 

The solution of equation (5), F(c) = 0 was found for the following cases; 

(a) Fixed resonant frequency of the diaphragm (q) for various <t> at 
the operating conditions of the fluidic generator. 

(b) Fixed <j> and varied r for these three operating conditions. 

The output of these results are shown in the following table. 
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It Is evident from the numerical results that 


(a) The natural Frequency of the system Is controlled by resonance 
conditions of the tube 1.e., wt/a * n/2, 3n/2 

(b) Oscillations are always damped 

(c) To keep the system In operation, forcing of the system Is 
needed. 
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V. DISCUSSION AND CONCLUSIONS 


The experimental work performed to develop the mathematical model to 
predict generator performance has shown that the mode of cavity operation Is 
governed by the standing acoustic resonance phenomenon in the tube. In this 
mode of operation, the annular jet Is periodically swallowed and subsequently 
discharged by the resonance tube at the fundamental tube resonance frequency 
(Reference 9). The basic mechanism of conversion of flow energy to acoustic 
power (and subsequently to mechanical motion) by the fluidic generator is 
critically dependent on the jet mass flow and the geometry of the gap between 
the nozzle exit and the resonating cavity. The basic understanding of the 
physical processes in this region under a simulated high altitude flow 
configuration is essential for designing the fluidic generator with optimum 
performance. More work is needed in this area to improve the efficiency of 
the conversion of fluid energy to electric power under high-altitude flight 
conditions. 

The comparison of laboratory test results with the flight test data 
(Figure 11) shows that the existing laboratory tests do not correctly simulate 
the performance of the fluidic generator at altitude. A close look at the 
flow field in the annular jet showed that the flow becomes supersonic inside 
the generator in flight, as compared with subsonic flow in the existing 
laboratory simulation facilities. To correctly simulate flight, a vacuum 
facility is needed to simulate both the high-altitude mass flow through the 
fluidic electric generator and the static flow pressures in and around the 
fluidic generator. 
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APPENDIX A 

CALCULATIONS OF LABORATORY CONDITIONS TO 
SIMULATE FLIGHT CONDITIONS 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


FORTRAN IV^PLUS V02-51 12{24J12 15-0CT-79 PAGE 1 

AfMi.FTN /TRZ BLQCKS/WR 



0001 

3002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 


0013 

0014 

0015 


0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 


0024 

0025 


C THIS PROGRAM DETERMINES CONDITIONS AT THE END OF THE 
C FLUIDIC GEN. GIVEN THE VELOCITY AND ALTITUDE OF THE MISSILE, 

C FROM THESE CONDITIONS AT THE END OF THE FLUIDIC GEN, 

C THE CONDITIONS UPSTREAM IN THE AIR STRAIGHTENING SECTION ARE DETERMINED 
C PS « STATIC SURFACE PRESSURE OF OGIVE AT PT. 2, 

CP - PRESSURE (N/M**2> AT DIFFERENT POINTS. 

C PO = STAGNATION PRESSURE (N/M*#2) AT DIFFERENT POINTS. 

C T « TEMPERATURE (KELVIN) AT DIFFERENT POINTS, 

C TO » TOTAL TEMPERATURE (KELVIN) AT DIFFERENT POINTS. 

CM * MACH NO. AT DIFFERENT POINTS. 

C PSI - F IN (PSD 

C TR = T IN <R) 

G PT. 1 ■ CONDITIONS UPSTREAM OF BQU SHOCK WAVE. 

C, 2 = CONDITIONS DOWNSTREAM OF BOW SHOCK (ACTUAL FLIGHT COND.) 

C JUST UPSTREAM OF FLUIDIC GEN. 

C PT, 3 * CONDITIONS JUST UPSTREAM OF FLUIDIC GEN. IN LAB. 

C APPARATUS, 

C PT. 4 » CONDITIONS IN THE AIR STRAIGHTENING SECTION OF LAB. 

C APPARATUS. 

DIMENSION HEAD(2) 

DOUBLE PRECISION P(4>»T(4),PSI(4)rTR(4>rP0(4)»T0(4)» 

1 CP(2)»CV(2)fM(4)»K»A3»A4» RHO(4)»PS * X 

DATA A3/1.266769AA4/182,41469/ 

DATA HEAD/4HAIR r 4HHEL♦/ 1 CP/3.5 f 2.5/»CV/2.5.1.5/ 

DATA XSAVE/1.0E+20/rPCQNV/l.45038E-04/ 

TYPE *»'INPUT VELQCIYT(M/S)r’ALTJTUDE(M)' 

ACCEPT *fVfZ 

CALL ATM76(Z»T(l)rP(l)?RHQ(1)) 

A * 20.046003 * SORT(T(1)) 

M(1) - V / A 

POM) ft P(l> * <1. + .2*M(1)#M(1))**3.5 
T0(1> « T(1) * <1, + .2*M(1)*M(1)) 

C 

C CALC. MACH NO. OF MISSILE. 

C 

A » 20.046803 * SQRT(T(i)> 

V = M(1) * A 
IF (M(l).GT.l.)G0 TO 5 
C 

C SUBSONIC FLIGHT 

C 

7(2) <= T( 1) 

V(2) <= P(l) 

M(2) m M( 1) 

P0(2) = P0(1> 

T0(2> b T0(1) 

PS = p(2) + .23 * (P0(2) - P(2)) 

GO TO 6 

5 CONTINUE 

C 

C CALC. OF CONDITIONS AT FLUIDIC GEN. FROM ATMOSPHERIC CONDITIONS 
C AND VELOCITY OF MISSILE.SAME PROBLEM AS CALC. DOWNSTREAM CONDITIONS 
r. OF A NORMA) .SHOCK WAVF. FROM .UPSTREAM CONDITIONS. 

C ' 

X * M(l) * M(l) 

P<2) b F'( 1) * <7.fX - 1.) / 6. 
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FORTRAN IV-PLUS V02-51 


PAGE 2 


0026 

0027 

0028 

0029 


0030 

003* 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 


0046 

0047 


0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 


/TR?BLOCKS/UR 

P0<2) s P<1> * <1.2*X)**3.5 * <P<1)/P<2))**2.5 

M<2) = SORT<(X+5.) / <7.#X~1.)> ' 

T<2) a T<1) * <7.*X~1.) * (X+5.) / 06.*X) 

T0<2) = T<2) * <1. + .2*M(2)*M<2)) 

CALC OF CONDITIONS JUST UPSTREAM OF FLUIDIC GEN. FOR 
THE LAB APPARATUS FROM THE ACTUAL UPSTREAM CONDITIONS 
REQUIRED CONSTRAINTS AREtDRIVING PRESSURE OF THE GEN. 

MUST BE THE SAME. 

I THE STAGNATION TEMP. OF THE ATMOSPHERE 
JUST BEFORE THE FLUIDIC GEN. <PT. 2) MUST EQUAL THE STAGNATION 
TEMP. AT PT. 3. 

DO 40 L«1r2 

K = CP(L) / CV(L) 

P<3> = 101325.0 

POO) a P0<2) + P<3) - PS 

M<3) = SORT< 2»/<K - 1.) * < <P0<3>/PO> >**< <K - l.)/K) - 1. 
STORE a ,5*<K - l.>*MO)*M<3) 

T0<3) a T0<2) 

T<3) a TOO) / <1. + STORE) 

EXPON a ,5 * <K + 1.) / <K — 1.> 

CRIT a A4/(A3#M<3)) * <1. + STORE)/<,5*<K+1 .> )**EXPQN 
X » Cl. / <.5*<K+1.)>)**EXPGN 
M<4) = X / CRIT 
P0(4) a POO) 

T0<4) = T0<3) 

P<4) = PO<4) / <1. 

T<4) a T0<4) / <1. 

PRINT OUT RESULTS. 


,5*<K~1. 
.5*<K-1< 


)*M(4)*M(4))**<K/<K~1)) 

)*M<4)*M<4)> 


30 

110 


120 

40 


TYPE 100.HEAD(L)»V»Z 

FORMAT</11Xr / GAS IN LAB APPARATUS? '»A4/2 Xf'VELOCITY = '.F8,2 

1 ' (M/S)'/2X»'ALTITUDE = '»F9.2*' <M)'»// 

2 31X»'STAG. STAG.'t/ 

3 lXr'PT. TEMP.<K) P <N/M**2) TEMP.(R) TEMP.(R) P (PSI)' 

4 ' P <PSI) MACH NO.') 

DO 30 1=1r4 

TR<I) = 1,8 * TCI) 

PSI(I) = PCONV * P<I) 

TRO a 1.8 * T0<I) 

PSIO = PCONV * PO(I) 

WRITE(A.110)I.T<I)fP(I)»TRO» TR(I).PSIO.PSI<I)»M(I) 

FORMAT <lXrI2r2X»F8.2»F13.2r2F10*2»2Fll*5»F12.6) 

X = <P0<2) - PS) * PCONV 
WRITE< 6»120)X 

F0RMAT<2X»'DRIVING PRESS. <P(4)-P<ATM.)> a ',F7.3) 

CONTINUE 

STOP 

END 
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PAGE 4 


FORTRAN IV-PLUS V02-51 
ATMl♦FTN /TR{BLOCKS/WR 

0001 SUBROUTINE ATM76<ZrTrPrRH0> 

C 

C SUBROUTINE ATM76 GENERATES THE ATMOSPHERIC CONDITIONS 
C TrPrRHO GIVEN THE ALTITUDE ACCORDING TO THE U.S. STANDARD 
C ATMOSPHEREr1976. 

C 


0002 


DOUBLE PRECISION TrPrRHOrMWTrGrRrROrZrH 

0003 


DATA R0/63S6766./r MWT/28.9655/rG/9.80655/r R/8314.32/ 

0004 


H * RO * Z / <R0 + Z) 

0005 


IF <H.LE.11000.)G0 TO 10 

0006 


IF <H.LE.20000.>60 TO 20 

0007 


IF (H.LE.32000.>GO TO 30 

0008 


TYPE *r'ALTITUDE IS TOO HIGH' 

0009 


RETURN 

0010 

10 

T = 288,15 - .0065 * H 

0011 


P “ 101325,0 * <288.15 / T>#*<<-l.>*G*MWT/<R*.0065>> 

0012 


GO TO 40 

0013 

20 

T » 216.65 

0014 


P « 22632.4 * EXP<<~1»)#G*MWT#(H-11000.>/<R*216.65) > 

0015 


GO TO 40 

0016 

30 

T = 216.65 + .001 * <H - 20000*> 

0017 


P = 5475.05 * <216.65 / T>**<G*MWT/R> 

0018 

40 

RHO = MWT*P / <R*T> 

0019 


RETURN 

0020 


END 


FORTRAN IV-PLUS V02-51 12:28:52 15-0CT-79 PAGE 5 

ATMl.FTN /TR:BLOCKS/WR 

PROGRAM SECTIONS 


NUMBER 

NAME 

SIZE 



ATTRIBUTES 




1 

♦C0DE1 000474 

158 


RWrIrCONrLCL 




2 

tPDATA 000136 

47 


RWrDrCONrLCL 




4 

♦VARS 

000050 

20 


RWrDrCONrLCL 




ENTRY POINTS 








NAME 

TYPE 

ADDRESS 

NAME 

TYPE 

ADDRESS NAME 

TYPE 

ADDRESS 

NAME 

ATM76 


1-000000 







VARIABLES 








NAME 

type 

ADDRESS 

NAME 

TYPE 

ADDRESS NAME 

TYPE 

ADDRESS 

NAME 

G 

R*R 

4-000010 

H 

R*8 

4-000040 MWT 

R*B 

4-000000 

P 

RHO 

R*B 

P-000010* 

RO 

R*8 

4-000030 T 

R*8 

F-000004* 

Z 
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I N?U? r 0tLOCIYT <H/S ) f ALT ITUDE < H) 
199.»4209. 


GAS IN LAB APPARATUS! AIR 
VELOCITY * 199.00 (M/S) 

ALTITUDE * 4209.00 <H> 


STAG. STAG. 


PT. 

TEMP.(K) 

P (N/M**2) 

TEMP.(R) 

TEMP.(R> 

P (PSD 

1 

240.81 

40000.78 

504.93 

469.46 

11.22968 

2 

240.81 

60000.78 

504.93 

469.46 

11.22968 

3 

270.73 

101325.00 

504.93 

487.31 

16.64199 

4 280.52 

DRIVING PRESS, 

114741.60 504.93 

> (P(4)-P(ATM.)> 1, 

504.93 

.946 

16.64199 


P (PSD 
8.70239 
8.70239 
14.69597 
16.64189 


GAS IN LAD APPARATUS! HEL. 
VELOCITY ■ 199.00 <H/S> 

ALTITUDE » 4209.00 (M) 


PT. 

TEMP.(K) 

P (N/M**2) 

STAG. 
TEMP.(R) 

TEMP.(R) 

STAG. 

P (PSD 

1 

260.81 

60000.78 

504.93 

469.46 

11.22968 

2 

260.81 

60000.78 

504.93 

469.46 

11.22968 

3 

266.91 

101325.00 

504.93 

480.43 

16.64199 

4 

280.52 

114741.62 

504.93 

504.93 

16.64199 

DRIVING PRESS, 
TTO -- STOP 

, (P(4)-P(ATM< 

*>> - 1, 

.946 



P (PSD 
8.70239 
8.70239 
14.69597 
16.64189 


>RUN ATH1 


INPUT VELOCIYT(M/S) pALTITUDE(H) 
176.9*4150.5 

GAS IN LAB APPARATUS: AIR 
VELOCITY * 176.90 (H/S) 

ALTITUDE = 4150.50 (M> 


STAG. STAG. 


PT. 

TEMP.(K) 

P (N/M**2) 

TEMP.(R) 

TEMP.(R) 

P (PSD 

1 

261.19 

60461.39 

498,17 

470.14 

10.73977 

2 

261.19 

60461.39 

498.17 

470.14 

10.73977 

3 

269.10 

101325.00 

498.17 

484.38 

16.21331 

4 276.76 

DRIVING PRESS. 

111786.15 498.17 

, (P(4)-P(ATM.)> - 1, 

498.17 

.517 

16.21331 


P (PSD 
8.76920 
8.76920 
14.69597 
16.21324 


GAS IN LAB APPARATUS: HEL. 
VELOCITY - 176.90 (M/S) 

ALTITUDE * 4150.50 (H) 


PT. 

TEMP.(K) 

P (N/M**2> 

STAG. 
TEMP.(R) 

TEMP.(R) 

STAG. 

P (PSD 

1 

261,19 

60461.39 

498.17 

470.14 

10.73977 

2 

261.19 

60461.39 

498.17 

470.14 

10.73977 

3 

266.10 

101325.00 

498.17 

478.97 

16.21331 

4 

276.76 

111786.16 

498.17 

498.17 

16.21331 


DRIVING PRESS. (P(4)-P(ATM.>) * 1.517 

TTO •— STOP 


P (PSD 
8.76920 
8.76920 
14.69597 
16.21324 


MACH NO. 
0.614676 
0.614676 
0.425257 
0.002850 


MACH NO. 
0.614676 
0.614676 
0.391152 
0.002585 


MACH NO. 
0.546016 
0.546016 
0.377305 
0.002548 


MACH NO. 
0.546016 
0.5A601A 
0.346783 
0.002315 
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INPUT VEL0CIYT(M/S) .ALTITUDE(M) 
165*7.3195.4 


GAS IN LAB APPARATUS! AIR 
VELOCITY » 185*70 (M/S) 

ALTITUDE = 3195.40 <H) 


PT. 

TEMP.(K) 

P (N/Mf*2) 

STAG. 
TEMP,<R> 

TEMP.(R) 

STAG. 

P (PSI) 

P (PSI) 

MACH NO. 

1 

267,39 

68397.19 

512,19 

461.30 

12.33323 

9.92019 

0,566492 

2 

267.39 

68397,19 

512.19 

481.30 

12.33323 

9.92019 

0.566492 

3 

275,04 

101325.00 

512.19 

495.06 

16.55401 

14.69597 

0.415938 

4 

284.55 

114135.07 

512.19 

512.19 

16.55401 

16.55392 

0.002792 

DRIVING PRESS, 

, (P(4)-P(ATM, 

,)) * 1, 

.858 





GAS IN LAB APPARATUS: HEL. 
VELOCITY * 185.70 <M/S> 

ALTITUDE ■ 3195.40 <M> 


STAG. 

PT. TEMP.<K> P (N/M**2) TEMP.(R) TEMP.(R) 

1 267.39 68397,19 512.19 .481.30 

2 267,39 68397.19 512.19 481.30 

3 271.32 101325.00 512.19 488.37 

4 284.55 114135.08 512.19 512.19 

DRIVING PRESS. <P(4)-P(ATM.>) » 1,858 

TTO — STOP 

>RUN ATM1 

INPUT VELOCIYT(M/S)rALTITUDE(M) 

203.4.1615,4 

GAS IN LAB APPARATUS: AIR 
VELOCITY -- 203.40 (M/S) 

ALTITUDE = 1615.40 (M) 

STAG, 

PT. TEMP.<K> P (N/M**2> TEMP.(R) TEMP.(R) 

1 277.65 83369.21 536.84 499.77 

2 277.65 83369.21 536.84 499.77 

3 284,45 101325.00 536.84 512.01 

4 298.24 119587.41 536.84 536.83 

DRIVING PRESS, <P< 4)-P(ATM.)) * 2.649 

GAS IN LAB APPARATUS? HEL. 

VELOCITY = 203.40 <M/S) 

ALTITUDE * 1615.40 (M) 

STAG, 

PT. TEMP.<K> P (N/M**2) TEMP.(R) TEMP.<R) 

1 277.65 83369.21 536.84 499.77 

2 277.65 83369.21 536.84 499.77 

3 279.11 101325.00 536.84 502.40 

4 298.24 119587,43 536.84 536.83 

DRIVING PRESS* (P(4)-P(ATM.>> = 2,649 

TTO — STOP 


STAG. 

P (PSI) 
15.53120 
15.53180 
17.34485 
17.34485 


P (PSI) 

12.09170 

12.09170 

14.69597 

17.34472 


MACH NO. 
0.608913 
0.608913 
0.453438 
0.002947 


STAG, 

P (PSI) 
15.53180 
15.53180 
17.34485 
17.34485 


P (PSD 

12.09170 

12.09170 

14.69597 

17.34472 


MACH NO. 
0.608913 
0.608913 
0.492382 
0.003261 


STAG. 

P (PSI) 
12.33323 
12.33323 
16.55401 
16.55401 


P (PSI) 
9.92019 
9.92019 
14.69597 
16.55392 


MACH NO. 
0.566492 
0.566492 
0.382521 
0,002533 
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CALCULATIONS OF MASS FLOW 
THROUGH THE FLUIDIC GENERATOR 


I. MASS FLOW CALCULATION THROUGH THE FLUIDIC GENERATOR 


It 1$ desirable to know the relationship between the absolute total 
pressure p 0 , the absolute cavity pressure p c , and the absolute skin pressure 
just outside the exit hole p 3 . To obtain these relationships, assume Ideal 
Isentroplc flow of air through a nozzle and multiply the mass flow rate 
expression by C v , the velocity coefficient: 

/*r Pipi <Pa/Po> 2/ V(P2/»o) Y " 1>/T ] 

m * C V A2\— . - . . — — 

Y_1 Ci-{p 2 /P a ) 2/Y (A 2 /Ai) 2 ] 

If state (0) Is the stagnation condition,then for A2 «Aj, the mass flow 
rate simplifies for air to 

m * C v A 2 J 7 Popo( P i)~ [i. ] 


Assuming the Ideal gas relationship P n » p Q RT n where R Is the gas constant, 


RTq W W 


m - (2.0546 CvApo/ZTolJ^J 7 [l- (^ 77 ] 


where A ~ In 2 , P 0 ~ Psla, T ~ °R, m * Ibm/sec. 

C v Is usually assumed a function of the Reynolds number of the orifice 

or nozzle. This Implies that the mass flow rate from a high pressure to a low 

pressure through some arbitrary small opening Is a function of the dimension¬ 
ed 

less quantities y * ratio of specific heats; R e = —; p/p 0 . The mass flow 

A |x 

rate Will also depend on the stagnation pressure and temperature, P 0 and T 0 . 
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The converging section of the bench setup allows for a smooth, gradual 
convergence of the air just before It flows Into the generator. The mass flow 
rate was calculated from the pressure ratio relative to this section, assuming 
that flow through this Section Is Isentroplc. This approximation should give 
an error In the range of ±5%, where the positive error Indicates that the mass 
flow rates calculated are greater than the actual mass flow rates. A simple 
procedure to get the actual flow rate Is to calibrate our setup using a 
flowmeter and obtaining the relationship of C v vs Reynolds number. The 
equation presently used to calculate the mass flow rate Is 



II. VELOCITY COEFFICIENT CALCULATIONS 

There are two velocity coefficients: C v j, the expansion from the ram 
pressure to the cavity, and C v n» the expansion from the cavity through the 
exit holes to the atmosphere. Each of these coefficients Is a measure of the 
effective area that Is reduced by the Influence of the boundary layer at the 
edge of each orifice. 

In the present Investigation, C v was studied as a function of Re only. 
The equation used to calculate C v j for the expansion from p 0 +p c 1$ 
m Vfo 

CvJ * 

Similarly, 


m YTo 
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The variables used to calculate C v j and C v h are described and their 


origins are noted In the table below. 



P 0 ,T 0 Stagnation conditions Experiment measurement 

m Mass flow rate Calculated from smooth 

converging section 


Aj,Ah Cross-sectional areas Drawings 

P s Skin pressure Experiment measurement 

P c Cavity pressure Experiment measurement 
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DERIVATION OF THE EQUATIONS OF 
FLUIDIC GENERATOR OPERATION 
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DERIVATION OF EQUATIONS OF FLUID GENERATOR OPERATION 


The operation of the fluidic generator has been divided Into five 
components consisting of the annular jet, the jet layer, the resonance tube, 
the cavity, and the holes (see Figures 1 and 7). A linear differential 
equation has been written for each of these five components that relate the 
mass flow and pressure fluctuations In different parts of the generator. 

Because of the linearity of the equation we look for solution of the form 
p' = p exp(a + 1w)t 

m' = m exp(cr + 1w)t (C.l) 

where w is the frequency of oscillations, cr the damping coefficient and p and 
m are the amplitude of the oscillations. A positive value of a indicates 
exponential growth at a particular frequency <o. From the homogeneous system 
of equations, the amplitude ratios can be obtained. These ratios contain a 
phase factor giving the phase differences between the various pressures and 
mass flow through the fluidic generator. 

The outcome of the system of these equations is a function a = a(w) which 
gives the damping or growth of acoustical waves of a given frequency and 
amplitude, and the phase relationship between pressures and mass flow at 
various locations. 

A. Annular Jet 

The annular jet is a convergent nozzle; the relation between mass flow 
and pressure ratio (total pressure before the contraction divided by the 
static pressure at the exit) is well known. In this analysis we have assumed 
that the exit pressure equals the cavity pressure, thus 
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Zy 2/y ( y -1)/y /2 \ 

ITIJ = C VJ Aj Po Po (P c /Po) Cl- (Pc/Po) 3 f° r (Pc/Po) ") 

My-1 y+1 


* 


C VJ Aj 


2y_ 

\ Y-l 


2 2/(y-1) y-1 
Po Po — 

Y+1 


~7 for Pc / P ° < (^) 


z \ y^t 


where p 0 , p c are the total and cavity pressures respectively. The second 
equation corresponds to shock conditions at the exit. The mass flow 
coefficient C v j incorporates possible viscous effects. Typical values of 
these coefficients are given in Table 1 to 4. 

Acoustical fluctuations in the cavity induce fluctuations of the mass 
flow which are given by the relationship. 


dmj 

dPc 


P’ 


dm. 


where ( )' denotes acoustical fluctuations and the derivative _^ is to be 

Wc 

evaluated at the operating point 

il! 

y+1 /Pc x Y 

1-— (-) Y/(Y-D 

2 V Po 7 pc ( 2 v 

- for - > (—) 

/ Pc \xdL Po V 1 

1 - \T) 1 


/dm 


m 


-)■ - 

dp c ' YPc 

W 


Po 


0 for 


pc 

Po 


<0 


y/(y-!) 
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The previous relation does not take into consideration the Inertial 
effects on the relationship between the acoustical pressure and mass flow. In 
order to Incorporate these effects, an Inertial term was incorporated into 
the equations. It Is assumed that the fluctuating pressure Is In part utilized 
to accelerate the mass flow through the nozzle. The resulting equation Is 



From dimensional considerations Ij will be a function of the operating 
conditions and the geometry of the jet nozzle. 


Ij * — fj (Mj, geometry) 

Aj 

where h Is the width of the ring jet, Aj its area and fj varies from a value 
of order one at Mj ** 0 to zero at Mj = 1. For the purpose of this analysis 
1 for Mj < 1 
0 for Mj > 1 


fj( M j) 


Seeking a solution as Indicated in equation (C.l), the equation for 
annular jet mass flow can be written as 


[Ij (o+iw) + 



Pc 


B. Pet Layer 

The jet layer is the region of interaction between the entering ram-jet 
and the lip of the resonant tube. In first approximation, the pressure 
inside and outside the ring jet layer can be considered uniform since the 
acoustical wavelength is much longer than the distance between the jet exit 
and the tube lip. The pressure inside (py) and outside (p c ) need not be the 
same instantaneously. Phase differences in their oscillations, however,will 
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result In a deflection force on the jet layer which will cause a radial 
displacement and therefore an oscillating mass flow Into the tube (my). 

Detailed analysis of this interaction is an extremely complex problem and 
is beyond the scope of current research. A phenomenological model was therefore 
developed in which the inertia and the curvature effects of the jet layer have 
been Incorporated. The present model assumes curved streamlines with constant 
radius of curvature as sketched. Furthermore, the jet width is assumed small 
compared with the radius of curvature of the layer and the ring jet. 


I 



Momentum conservation requires that the pressure differential be balanced 
by the inertia and the centrifugal force. The inertia term is given by 



where 



we find 





*2 d(l/Rj) 


= phWq — -- 
2 dt 


o 

where h is the jet width, W its perimeter, and q the jet exit velocity assumed 
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constant and Independent of x. The centrifugal term Is given by 



and finally the pressure force Is (Py - P c ) AW. 


Thus 


4 .) 


x d W , 

P h q-- + p _ h , ( p T . p c) 

2 dt Rj 

The radius of curvature Is related with the tube mass flow m' by the equation 

m* J 

2 T 

1/ R J - - - — 

A 2 pqW 

The relation between the tube mass flow, the pressure in the cavity and the 
tube is 

I # 

H d lily 2 (m) w 

— - + —To 55 P' - P' 

AW dt pA W z T c T 

where (m) w is the mass flow through the system at the operating conditions. 
For sinusoidal solution we find: 


r h 2(m) w -.a 

I “ (o+io)) + —-- my = P c - P T 

l aw p a 2 w 2 j 

C. Resonant Tube 

The most important aspect of the operation of the fluidic generator is 
the energy transfer in the resonant tube. An adequate model of the system has 
to incorporate its two functions as a) the resonant tube system and b) the 
"acoustical transmission line" phenomenon through which the energy is 
transferred to the diaphragm. Depending upon the acoustical wavelength, a 
variety of tube resonant modes can be excited. At low frequencies, i.e, the 
wavelength is much longer than the diameter of the tube, only the longitudinal 
modes are excited. These modes are characterized by plane waves moving back 
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\ 

* 

and forth through the tube. The difference In energy carried by the waves 
travelling in one direction compared to those moving in the opposite direction 
is the amount of energy transferred through the tube. The resonant 
characteristics of the tube are determined by the matching conditions at both 
ends which in turn determine the amount of energy transmitted. 

In the present model, the relationship between the pressure and the mass 
flow at the entrance of the tube Is a function of the amount of energy 
transferred through the tube. 



We consider the superposition of right and left propagating acoustical waves. 
Their amplitude and phase at the diaphragm end are related by the reflexion 
coefficient R(w) which in general is a complex function of the frequency. 

In general, if the pressure field associated with the right propagating 
wave at x=o is 
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then at any x 


P R » P„(t) 


PR(t.*) ■ PR (t- —) 
a 

where a Is the speed of sound. Similarly for the left propagating wave 

p L (t,x) = PL(t+-) 
a 

At the diaphragm end, the reflexion coefficient R Is defined as 

R . p L(t + l/a) 

P R (t- L/a) 

For wave forms given by 

p R (t) = Pr exp[(a+1w)t] 

P L (t) = P L exp[(cr+iw)t] 

we find 

p L 

R =— exp [2(o+iw) L/a] 

P R 

and consequently at x=o the pressure fluctuation is given by 

p T - ( P R + Pl) exp[(<j+1w)t] = p R exp[(o+i w )t] [1+R exp {-2(a+iw)L/ao}] 
Associated with the right and left propagating waves there is an acoustical 
velocity field given by 

Ur = P R /pa 
U|_ = - p L/pa 

Therefore, the velocity at the entrance of the tube is given by 
p R 

U T = —- exp[(crti{ji))t] {1-R exp(-2(o+ia)) L/a)} 

Pa 

The mass flow mj = p Uy Ay is 

mj = my exp[(a+iw)t] 

Pr 

ihy - — Ay {1-R exp [-2(a+iw)L/a]} 
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The relationship between mj and pj Is 

m T Aj 1 . R exp C-2(a+1w) L/a] 

PT a 1 + R exp J>2(o+1w) L/a] 

The reflection coefficient R Is determined by the Interaction of the 
acoustical field with the diaphragm* The diaphragm can be modelled as a 
second order mechanical system. If xo Is its position then Its notion is 

d 2 x D dx D 

described by the equation P'p V m D — + r D —- + k 0 x D 

dt2 dt 



where mg and kg are the mass and spring constant of the diaphragm-reed 

combination, and fp is the friction coefficient which also incorporates the 

effects of the energy transmitted into the electrical part of the system. The 

dxo 

fluid velocity at the diaphragm is —which results In a velocity at the end 

dt 


of the tube IT, 


Ad Ad dxo 

U f = — U' -- 

£ Ay d At dt 


Since p'j: at the end of the tube equals p^, we find 


The tube equation is given by 
P E 


> a t /a d 


d 2 xo dXD 

mg - + fo — + k D x D 

dt dt. 


dx D 
A D —- 

dt 


U E 


pa 


HR 

HR 
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The reflection coefficient can be written as: 


paAp dxp d 2 xp 

Ayfp dt dt 2 


paA d dxo d 2 Xp 

(1+ -) -+ mo—-- 

Ajfp dt dt 2 


For sinusoidal oscillations, one can write 

.2 

paAp 

( 1 - -) i + mD z 2 * ko 

ATfD 


paAp 

( 1+ -) z + mp z 2 + kp 

A T fD 


where z = c+1u> 

Introducing the nondlmenslonal parameters 

paAp 

ATfp 

_ m[) a 
11 ’ fD L 

b ,JE I 

\ m a 


The term coefficient can be written as 


R = 


Lz p Q 2 -i 

W+n— rit(--) ] 

a L Lz/a J 

Lz r D 2 ^ 

1+4>+t) i+(-) j 

a L Lz/a J 


+ kpxp 


+ kpxp 
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0. Cavity 

Both the jet and the tube are located inside a cavity which plays a role 
In the acoustical properties of the system. For the frequencies of fluidic 
generator operation, this cavity is operating as an accumulator because the 
wavelength Is very large compared with the characteristic dimension of the 
cavity. The acoustical pressure Inside the cavity can also be considered 
uniform. Furthermore, the characteristic time of the oscillations is short 
enough to neglect the heat transfer from the fluid to the wall. Therefore the 
cavity can be modeled as an isentroplc accumulator for which the relation 
between the mass flow and cavity pressure Is given by 
dp c yRT c , 

dT ~ 

where y is the ratio of specific heats, R the gas constant, T c the gas 
temperature and V c the cavity volume. 

This equation gives for oscillatory waves 

- yRT c f 
(a+1w) p c = -— m c 
Vc 
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E. Holes 


The equation modeling the behavior of the exit holes is Identical to that 
governing the behavior of the annular jet. It should be noted that for the 
exit holes, the acoustical fluctuations depend on the total pressure upstream 
of the orifice rather than at the exit pressure; thus: 

I v-1 y/y- 1 

VcJj[.« '],2»(1) 

V Y-l 'Pc' L 'Pc' J Pc Vl' 

I Y/Y-l 

2 Ps / 2 \ 

* Cvh Ah Jy(—) p c Pc ; — <1-1 

V y+1 Pc V +1 ' 

Therefore the linearized equation Is 
d ifiH« /dp C \ 


IH 


H‘ / dPc \ 
Vdinu / 


dt VdifiH' H 


with 


/ dPc\ / dl ^H\ 


( tTi(Y-l)Pc 

2-(Ps/Pc) 
l-CPs/Pc)* -1 

; Ps/Pc> | 

(-) 

2 Y 

/ 9 \Y 

1 * L 

V Pc 


; p s /Pc< 1 

© 


I H = — f H (M ex , geometry); f« (M ex , geometry) = 1 
ndn 

f H has been assumed independent of the exit Mach number. 
Thus for sinusoidal oscillations, we find 
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COMPUTER PROGRAM FOR MODELING THE 
FLUIDIC POWER GENERATOR 
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Table 1. Data for All Exit Holes, Open Area 
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Table 3. Data for 3/4 Exit Holes, Open Area 
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Figure 1. Fluidic Electric Power Generator 
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Figure 2. Experimental Setup for Fluidic Generator Bench Test 
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Figure 3. Instrumentation of Fluidic Electric Power Generator 
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Figure 7. Calculation of Bench Flow Conditions 
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Figure 8. Comparison of Labora ry Measurements with Flight Data 
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Figure 11. Comparison of Laboratory Measurements with Flight Data 




Figure 12. Mass Flow as a Function of P./P. for Various Exit Areas 
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Figure lb. Velocity Coefficient vs Reynolds Number for Various Exit Areas 
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Transducer 1 Resonance Tube 
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Autocorrelation 
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Figure 16. Pressure Transducer Output 
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Figure 18. Oscilloscope Traces of Fluidic Generator 
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Figure 19. Transient Response of Diaphragm-Reed Subsystem 
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Figure 20. Modeling of Fluidic Electric Generator 







